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Abstract. Cellular automata are widely used to model real-world dynamics. We show using the Domany- 
Kinzel probabilistic cellular automata that alternating two supercritical dynamics can result in subcritical 
dynamics in which the population dies out. The analysis of the original and reduced models reveals gen- 
erality of this paradoxical behavior, which suggests that autonomous or man-made periodic or random 
environmental changes can cause extinction in otherwise safe population dynamics. Our model also real- 
izes another scenario for the Parrondo's paradox to occur, namely, spatial extensions. 

PACS. 02.50.Ga Markov processes - 05.50.+q Lattice theory and statistics (Ising, Potts, etc.) - 87.23.Cc 
Population dynamics and ecological pattern formation 

1 Introduction the Domany-Kinzel (DK) model, which is a two parameter 
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d ' family of Markov processes on a one-dimensional lattice 

_ _ i 

Ecological and sociological dynamics are often described with discrete time H31- In this paper, we report a coun- 
by systems of locally interacting agents. Cellular automata terintuitive phenomenon of the DK model: particles even- 
are broadly used for modeling such dynamics to charac- tuall y die out when tw0 supercritical DK dynamics alter- 
terize, for example, survival probability, percolation, and natc with some appropriate orders. This behavior is robust 
critical phenomena, which are relevant to real situations against parameter changes. We also analyze the reduced 
P. Among the class of probabilistic cellular automata is dynamics such as the pair approximation and a canonical 
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model to guarantee that this phenomenon is preserved in 
much simpler models. As a generalization, dynamic en- 
vironmental changes can extinguish a population even if 
the snapshot dynamics is supercritical at any given mo- 
ment. These alternating DK dynamics also realize a new 
scenario for the Parrondo's paradox to occur, that 

is, introduction of the space. 

2 DK model 

In the DK model HUH], each site either accompanies a 
particle (denoted by •) or is empty (denoted by o) at 
any instant. The space can be identified with a subset 
of the set of integers Z, and let („ c Z be the set of 
the sites that have particles at discrete time n € Z + 
= {0,1,2,...}. The stochastic evolution rule at each site 
x £ Z is independently described by P(x £ £n+i|£n) = 
f{\tn n {x- l,x+l}\) where /(0) = 0, /(l) = Pl , /(2) = 
P2, and (pi,p2) £ [0, l] 2 . In other words, the probabil- 
ity that a particle emerges is determined by the number 
of the particles in the nearest neighborhood in the pre- 
vious time, as shown in Fig. ^ Each realization of the 
spatiotcmporal process is expressed in the form of a con- 
figuration £ e {0, 1} S = X with S = {s = (x,n) € 
Z x Z + : x + n = even}. The region of the supercriti- 
cal parameter sets (p\,P2) for which particles survive for 
infinite time with positive probability can be numerically 
obtained, and it occupies an upper-right area in the p±-p2 
space |3El[7j. The DK model is equivalent to the directed 
bond percolation model on a square lattice when (pi,P2) 
= (p, 2p — p 2 ) and to the directed site percolation model 
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when pi = P2 = P lUdllZ Another special case is Wol- 
fram's rule 90 deterministic cellular automaton [Q which 
is realized with (pi,P2) — (1 ; 0). The simplicity of the 
DK model enables us to investigate interesting proper- 
ties from the viewpoint of statistical physics and applica- 
tions, such as quasistationary particle density [51IU1I10II11I 
112) . critical phenomena and phase transitions 23,8; K?lllOI 
ITT1|12II13II14| . survival probabilities 03, and duality [Tol 
EttH 

Let us denote by P n (-) the probability that an event 
occurs at time n. Here an event means a state of consec- 
utive sites, or a sequence of • and o. For clarity, we often 
plot trajectories in the two-dimensional space spanned by 
the order parameters defined with a-iin) = P n (**) and 
ai{n) = P„(»o) + P„(o«). With a Q (n) = P„(oo), it fol- 
lows that a\(n) > 0, 0,2(71) > 0, and ai(n) + 02(11) = 
1 — ao(n) < 1. The origin (aj., 0,2) = (0, 0) is an absorbing 
fixed point corresponding to the population death. In the 
following numerical simulations, the lattice size is 10000, 
and the periodic boundary conditions are assumed. 

With some initial conditions, trajectories of the DK 
model are shown in Fig.^a) for (pi,P2) — (0.52, 1) (thin 
lines) and for (0.76, 0.76) (thick lines). The DK dynam- 
ics corresponding to these parameter sets are termed dy- 
namics A and dynamics B, respectively. When P2 = 1, 
particles emerge or die only at kinks where • and o face 
each other. In this case, the dynamics of kinks are iden- 
tical to the coalescing random walk, and the entire space 
is eventually occupied by particles with a positive prob- 
ability if and only if p\ > 0.5 0. Therefore, dynamics A 
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is supercritical. On the other hand, the DK model with 
Pi = P2 is equivalent to the directed site percolation. Re- 
stricted onto this line, p\ = P2 = 0.75 is a mathemati- 
cally rigorous upper bound for the subcritical regime 
whereas the critical value is numerically estimated to be 
about pi — p2 = 0.7055 |2(JU21| . Because of the attractive- 
ness (pi < P2), the natural intuition that more particles 
are likely to survive with larger p\ and P2 actually holds 
|7lll5j . Therefore, dynamics B is also supercritical. Ac- 
cordingly, trajectories of dynamics A converge to the all 
• state, and those of B converge to the stochastic stable 
fixed point (01,02) = (0.39,0.42). 

3 Population Death in Alternating DK 
Dynamics 

Next, we alternatively apply A and B. A typical trajectory 
is shown in Fig-EJb) with the Bernoulli initial distribution 
with density 0.5, which yields (oi(0), o 2 (0)) = (0.5,0.25). 
Surprisingly, particles eventually die out. This behavior is 
not sensitive to the choice of initial conditions. It also per- 
sists against changes in p\ or P2 as far as the individual 
dynamics are not extremely supercritical and the stable 
stochastic fixed points for the two systems are separated 
enough. Especially, extensive numerical simulations sug- 
gest that this population death is enhanced when one of 
the component dynamics is nonattractive, or p\ > P2- 

An important cause for the population death is how 
the trajectories of dynamics A and those of dynamics B 
cross. As shown in Fig.[2Ia), if a state in the 01-02 space 
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evolves along a trajectory of dynamics A, in terms of dy- 
namics B, the state gradually slides down to trajectories 
associated with initial conditions with fewer particles. In 
other words, from the viewpoint of dynamics A (resp. B), 
the population once decreases under dynamics B (resp. 
A) before it revives and reaches the nontrivial stable fixed 
point. Therefore, by switching the dynamics between A 
to B before the population effectively starts to grow, the 
number of the particles gradually decrease to zero. Sur- 
vival results if A or B is applied long enough before switch- 
ing to the other. To demonstrate this, we confine ourselves 
to the cases in which a block of k A's and k B's are al- 
ternatively applied, which we denote by A k B k . As shown 
in Fig.^a), the population is more likely to survive as k 
increases. 

The population death by alternation is an example of 
the Parrondo's game in which a combination of two losing 
(winning) stochastic games can counterintuitively end up 
with a winning (losing) game 4 5 6, 22 . In this context, 
the results in Fig. 01a) agree with those for the original 
Parrondo's game in which the paradoxical effect becomes 
small as k is raised [B]. A more general concern is how 
the arrangement of A and B affects the upshot. Since it 
appears quite difficult to derive the optimal ordering of 
A and B among all the possible sequences pOES), we only 
deal with some representative cases. 

The population dynamics when a chain of A is periodi- 
cally punctuated by just one B, which is denoted by A k B, 
are shown in Fig. EHb). This figure together with addi- 
tional numerical simulations suggests that the paradoxical 
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effect is most manifested, or the population dies out most 
rapidly, with k = 4. This is presumably because dynam- 
ics B correspond to the critical line of the attractiveness 
(pi = P2)- For this reason, in an upper-left region in the 
0,1-0,2 space, an application of the near-nonattractive B 
kills more particles when there exist more particles. This 
view is supported by Figs.yfb) and|2[c) in which we com- 
pare the dynamics with AB and those with A 4 B. Then, 
the convergence to (01,02) = (0,0) is accelerated by a 
larger A: in a small k regime. However, with a much larger 
fc, the number of particles changes little for most of the 
time (Fig. |2Jd)). In this regime, the population death is 
slowed down as k increases. 

For sequences in the form AB k , the parity effect is 
manifested. As shown in Fig. Etc) , the population death 
is faster when k is even. This is again because dynamics 
B is nearly nonattractive. As is prominent in nonattrac- 
tive DK dynamics, the motion in the 01-02 space under 
dynamics B is somewhat sensitive to the current state. 
More specifically, simple repetition of B yields a damped 
oscillation in the early stage. Therefore, if the initial state 
is located in a upper-left region, the number of particles 
drops more when B is repeated even times before being 
interrupted by one A. 

The random arrangement of A and B is also of in- 
terest 1 1IGI22| because real environments can be random 
rather than perfectly periodic. To mimic simple random 
environments, we choose A and B independently at each 
time step with probability r and 1 — r, respectively. Ob- 
viously, the population death does not occur with r = 
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or with r = 1, which prescribes the sequence purely of 
B and that of A, respectively. Figure 0] and the exten- 
sive parameter search reveal that the paradoxical effect is 
maximized when r = 0.2. This value coincides with the 
optimal mixing ratio for the family of deterministic se- 
quences investigated above, namely, A 4 B. 

It is also essential for the paradox that population 
change rates are proportional to the population size as 
shown in Fig. [3] (crosses). Owing to this property, the size 
of the population exponentially shrinks to a very small 
level (Figs. |3 El El- Then, particles become extinct in fi- 
nite time because of stochasticity and the absolute stabil- 
ity of the fixed point (aj., 02) = (0, 0). If change rates are 
too high even for minute population mass, a trajectory 
that happens to have approached the origin more likely 
escapes the vicinity of the origin to avoid the population 
death. 

4 Pair Approximation 

To take a closer look at the paradox, we analyze the de- 
terministic dynamics derived by the pair approximation 
of the DK model, which we call the PA dynamics 
I11I12| . In the pair approximation, any events at two sites 
separated by a distance more than one are supposed to be 
independent of each other. For example, P„(»|o«) = P„(»o 
• )/P„(o«) is approximated by P„(»|o) = P„(»o)/P n (o), 
where P n (-|-) denotes the conditional probability. Accord- 
ingly, probabilities of any events involving three or more 
consecutive sites are decomposed into one- or two- site 
probabilities. With this approximation, the two-dimensional 
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PA dynamics are written as follows: which has fixed points b\ = and b\ = (2p\ — l)/(2p\ — p 2 ) 

[5lll0lllllll2| . Let us pick two mean-field dynamical sys- 



a 2 {n + 1) = P 2 Pn(» ••) + PlP2{Pn(* • o) + P n {° • •)} 
+ Pi{Pn{» o •) + P„(o • o)} 



{2p 2 a 2 (n) +p 1 a 1 {n)) 2 p?ai(n) 2 



terns so that their nontrivial fixed points are positive and 
stable, with 2p\ > p 2 and p% > 1/2 satisfied. Then, when 



(1) 



46i(n) 4&o(^) ' two mean-field dynamics alternate, the particle density 

la 2 (nf 
h(n) 



2a 2 (n) 2 

a\(n + 1) = p 2 {l — p 2 ) bi(n) just moves between these two fixed points in the 



i o \ a i( n )°3( n ) lone; run. Accordingly, the population never dies out, and 

+ (Pi +P2 - 2pip 2 ) — — 6 6J ' 

bi(n) 

i \2 / \ / \ uo paradoxical phenomenon occurs. 

2b {n)bi(n) b (n) 

where bi(n) = a 2 (n) + ai(n)/2, bo(n) = 1 — b\{n) and 

/ \ , , s ( \ rp_ • . r i-i J3 * j 5 Canonical Model 

ao{n) = 1 — ai(rij — 02 (n). Irajectories ot the PA dy- 
namics are shown in Fig. Ela) for two sets of supercritical 

/ \ rr, k« ,i /, n . ,. \ , To generalize the Parrondo's paradox found for the DK 

parameter sets: {pi 7 p 2 ) = (0.52, 1) (thm lines) and (0.66, 

r, i i- \ T i ■ , i T-i- RT \ , i and PA dynamics, we construct a simple canonical model 

0.66) (thick lines). In accordance with lig. |2Ja), the m- 

,. . , , , . , , , „ , . . x with dimension two, which is the presumed minimal de- 

dividual PA dynamics own stable fixed points near (0, 1) 

i irs ono r> 1 a r\ tt i . tt n*vi \ n gree of freedom for the paradox. As we have mentioned, 

and (0.308, 0.145). However, as shown m Fig. |p[b), the 

. ; . ,. . ,. ,, A1J , , the relevant features of the DK and PA dynamics can be 

population dies out wncn tney alternate. Although the 



supercritical parameter region of the PA dynamics devi- 
ates from that of the DK counterparts, the results for the 
PA dynamics qualitatively agree with those for the DK 
dynamics shown in Fig. [21 

The Parrondo's paradox is unlikely to happen in one- 
dimensional systems since they lack auxiliary dimensions 
that counteract the seeming tendency of population in- 
crease. To demonstrate this, let us imagine the simplistic 
mean-field approximation in which a joint probability is 
approximated by a product of single-site probabilities (e.g. 
P„(o») = P„(o)P„(»)). The approximate one-dimensional 
system is written as 



summarized as follows. 



(i) Trajectories of dynamics A and those of dynamics B 
transverse in the way as shown in Figs. Ufa) and||Ja). 
More specifically, in the a\-a 2 space, the slope of a 
trajectory of dynamics A (thin lines) is less negative 
than that of a trajectory of dynamics B (thick lines) 
at the crossing point, at least in a certain region. 

(ii) Each of A and B is not applied too many times suc- 
cessively. In other words, k in the sequence A k B k , 
A k B, or AB k should be small enough, as explained 
with Figs. Eld) and^a). 

(iii) Population change rates are proportional to the pop- 
&i(n + l) =p 2 b 1 {n) 2 + 2p 1 b 1 (n)(l-b 1 (n)), (3) ulation size. To weaken the condition may result in the 
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same conclusion just with a different convergence rate. 
Here we assume this linearity for our canonical model. 

The paradox also relies on the following implicit as- 
sumptions. 

(iv) Dynamics A and dynamics B have sufficiently sepa- 
rated nontrivial fixed points. 

(v) The origin is the deterministically unstable but stochas- 
tically reachable fixed point for both A and B. 

Figure Efa) and EJa) further indicate that the origin 
and the nontrivial fixed point are connected by applying 
B infinitely many times (thick lines). However, it is not 
true for A (thin lines) because A is nongeneric in the sense 
that all the points on the ai axis are fixed points. Actually, 
no point with a\ — and < a<i < 1 is realizable because 
it would mean that two consecutive sites take state •• and 
oo with positive probabilities but not »o or o». The total 
size of the boundaries between clusters of • and those of 
o, or a\(n), determines the population change rates |7]. It 
declines to zero as a point in the a\-a2 space approaches 
the ai axis. 

In fact, we have chosen (pi,p2) = (0.52, 1) for dynam- 
ics A just because the obtained DK dynamics is rigor- 
ously supercritical. The paradoxical dynamics appear ro- 
bustly against changes in pi and P2, which can make the 
dynamics generic. With this in mind, we construct a two- 
dimensional continuous-time system that satisfies the con- 
ditions listed above. We propose to alternate two dynam- 
ical systems: 
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dynamics A 

V = Ay(l - y), 

and 

dynamics B 

x = Xx(l — x), 

is = -y, 

where < x, y < 1. The properties (iii), (iv), and (v) 
are obviously satisfied, with (iii) also supported by Fig. [5] 
(squares). Both dynamics A and B have a saddle at the 
origin. The point (0, 1) of A and (1, 0) of B are stable equi- 
libria, and each of them is connected to the origin by a 
heteroclinic orbit. The property (i) is satisfied if < A < f . 
To guarantee (ii), we set A = 0.3 and the duration of each 
dynamics equal to 0.15. Figure {7\ summarizes flows of the 
individual dynamics (thin lines) and those of the alternat- 
ing dynamics (thick lines). We again observe the paradox 
that the alternating dynamics lead the state toward the 
origin. 

6 Conclusions 

We have shown using the DK model and its simplifications 
that mixtures of two supercritical dynamics can yield sub- 
critical dynamics in which the population dies out. This 
counterintuitive behavior occurs if individual component 
dynamics have at least dimension two and satisfy certain 
criteria. The property (i) is characteristic of the DK or 
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the canonical model, and it agrees with some natural oc- which should be compared to PCA with m = 3. One 

casions but not with others |23j . The other four require- reason is that outcomes depend not only on the number 

ments do not seem to spoil the reality. The properties (iii) but also on the arrangements of particles in a neighbor- 

and (v) are satisfied when production rates are primarily hood [Sj. For instance, it is easy to verify P n (*\ • ° # ) 7^ 

proportional to the population mass, which is quite com- P n ( u \ • ,0 )- More importantly, the alternating DK dynam- 

mon for ecological and social systems [23- Periodical and ics are not special cases of finite-range PCA. To illustrate 

random environmental changes comply with (ii) and (iv). this, let us consider ABAB . . .. In Fig. the state of site 

Such changes may be also caused by continual, periodic, a, which we write £(a) depends on £(/), and 

or random human control of a system with the aim of while £,{h), and put together determine £(&). In 

moving the stable fixed point to more desirable one. How- 3-neighborhood PCA, there exists no intermediate layer 

ever, our results indicate that environmental changes or of sites such as c, d, and e. Therefore, once £(/), £(ff), 

oddly managed control measures can cause a total disaster £(/i), an d are given, £(a) and £(6) are independent, 

even if the system instantaneously stays in a supercritical On the other hand, in our model, £(a) and £(&) are par- 

'good' regime all the time. The other way round, there is tially correlated, or correlated even conditioned by £(/), 

a general expect that a situation that is subcritical at any £(<?), £,{h), and £(«). This is because both £(a) and £(6) 

moment can be changed into a supercritical one with ap- depend on £(<i). By the same token, the infinite-range cor- 

propriate controls, which is originally illuminated by the relation is generated just after single application of AB, 

Parrondo's paradox [4"U5UuTl22| . In the context of the Par- which prohibits use of powerful duality equations [1611171 

rondo's paradox, our model provides another mechanism I18j . In this sense, our model stiupulates a class of infinite 

of its occurrence in addition to inhomogeneous game rules particle systems different from ordinary PCA. However, 

or players with memory 5 , namely, spatial extension. on the analogy of the Parrondo's paradox, the phenom- 
ena reported in this paper may hold for PCA and more 

Lastly, we can regard a block of sequence of A and 

general alternating dynamics with general neighborhood 

B, such as A k B and AB k , as a transformation done in 

sizes. 

just one step. By doing so, the alternating DK model 
seems similar to m-neighborhood probabilistic cellular au- 
tomata (PCA) with m > 3. For PCA, phase diagrams have 

been studied in simple cases where the dynamical rule de- We thank K. Sato for his helpful comments. This study is sup- 

pends only on the number of particles in the neighborhood ported by the Grant-in- Aid for Scientific Research ( JSPS Fel- 

with m = 3 |13U14j . However, the model proposed here is lows) and the Grant-in- Aid for Scientific Research (B) (No. 12440024) 

more complex even with the simplest sequence ABAB . . ., of Japan Society of the Promotion of Science. 
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Figure captions and PA dynamics and is denned to be x 2 + y 2 for the 

canonical dynamics. 

Figure 1: Schematic diagram showing the DK proba- 
bilistic cellular automaton. Figure6: (a) Trajectories of the PA dynamics for { Pl ,p 2 ) 

(0.52,1) (thin lines) and those for {pi,p 2 ) = (0.66,0.66) 

Figure 2: (a) Trajectories of the DK model for dynam- ( thick lincs ). ( b ) Population dynamics in the alternating 
ics A with ( Pl ,p 2 ) = (0.52, 1) (thin lines) and those for pA dynamics starting from ( &1 (0), a 2 (0)) = (0.5,0.25). 
dynamics B with (pi,p 2 ) = (0.76,0.76) (thick lincs). The 

Figure 7: Trajectories of dynamics A and those of dy- 

other panels show population dynamics when we repeat 

„ . , , . ,„ , , „ ,, n m , , , „ namics B of the canonical model (thin lines), superim- 

(b) AB, (c) A 4 B, and (d) A 30 B. The initial conditions for v h 1 

posed by those of the alternating dynamics starting from 

(b, c, d) are (oi(0), a 2 (0)) = (0.5,0.25). 

(x,y) = (0.5,0.5) (thick line). 
Figure 3: Dynamics of the population size when we Figure 8: Alternating dynamics with sequence AB corn- 
repeat (a) A k B k with k = 1 (thinnest line), 2, 3, and 4 pared with standard 3-neighborhood PCA. 
(thickest line), (b) A k B with k = 1 (thinnest), 2, 4, 15, 
30 (thickest), and (c) AB k with k = 1 (thinnest), 2, 3 
and 4 (thickest). In (b), the lowermost line corresponds to 
k = 4. In (c), the upper lines, which are nearly superim- 
posed, correspond to k = 1 and 3, whereas the lower lines 
correspond to k — 2 and 4. 

Figure 4: Population dynamics when A and B ran- 
domly appear with probability r and 1 — r, respectively. 

(a) r = 0.02, 0.05, 0.1, 0.2 (from upper to lower lines), and 

(b) r = 0.2, 0.4, 0.5, 0.75 (from lower to upper lines). 

Figure 5: Population change rates in terms of the pop- 
ulation size for the DK (crosses) , PA (circles) , and canon- 
ical (squares) dynamics. The change rates for the DK and 
PA dynamics are measured by the Euclidean distances of 
two points with unit time difference in the a\-a 2 space. 
The population size is equal to a\(t)/2 + a 2 (t) for the DK 
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